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Abstract 

We present preliminary measurements of the CP asymmetry parameters and CP content in B° — > 
K + K~K® decays, with cpK® events excluded. In a sample of 227 M BB pairs collected by the BABAR 
detector at the PEP-II B Factory at SLAC, we find the CP parameters to be 5 = -0.42±0.17±0.04 
and C = 0.10±0.14±0.06, where the first error is statistical and the second is systematic. Extracting 
the fraction of CP-even final states from angular moments f even = 0.89 ± 0.08 ± 0.06, and setting 
C = 0, we determine sin 2(3 = 0.55 ± 0.22 ± 0.04 ±0.11, where the last error is due to uncertainty 
on the CP content. 
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1 INTRODUCTION 



In the Standard Model (SM) of particle physics, the decays B° —> K + K~K® [1] are dominated 
by b — > sss gluonic penguin diagrams [2]. CP violation in such decays arises from the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing mechanism [3]. 

The time-dependent CP asymmetry is obtained by measuring the proper time difference At 
between a fully reconstructed neutral B meson (Bqp) decaying into K + K~ K®, and the partially 
reconstructed recoil B meson (-B^ag)- The asymmetry in the decay rate f+(f-) when the tagging 
meson is a B° (B°) is given as 

-|At|/r s0 

f±(At) = — [ l±Ssin(Am d At)=F C cos (Am d At)], (1) 

4r B o 

where r^o is the B° lifetime and Arrid is the B°-B° mixing frequency. The parameters C and S 
describe the magnitude of CP violation in the decay and in the interference between decay and 
mixing, respectively. In the SM, we expect C = because there is only one decay mechanism and 
direct CP violation requires amplitudes with different CP-violating phases. If K + K~K® decays 
proceed through a P(S)-wave leading to a CP-odd (even) final state, we expect S = (— ) sin 2/3, 
where sin 2/3 = 0.731 ± 0.056 [4, 5]. 

However, contributions from physics beyond the SM could invalidate these predictions [6]. Since 
b — > sss decays involve one-loop transitions, they are especially sensitive to such contributions. 
Recent results in decays of neutral B mesons through the (f>Kg intermediate state are inconclusive 
due to large statistical errors [7, 8]. 

A more accurate CP measurement can be made using all decays to K + K~K® that do not 
contain a (f> meson. This sample is several times larger than the sample of 4>K® [9, 10, 11], but 
the CP content of the final state is not known. In this paper we present measurements of the CP 
content in K + K~K® decays using an angular- moment analysis and cross-check it with an isospin 
analysis [10]. We update our measurement of the CP asymmetry parameters and extract the SM 
parameter sin 2/3 with almost twice the statistics of the previous BABAR result [11]. 

2 EVENT SELECTION 

This analysis is based on about 227 million BB pairs collected with the BABAR detector [12] at the 
PEP-II asymmetric-energy e + e~ storage rings at SLAC, operating on the T(4S) resonance. 

We reconstruct B mesons from K® — > ir + ir~ and K ± candidates. Charged kaons are distin- 
guished from pions and protons using energy-loss (dE/dx) information in the tracking system and 
from the Cherenkov angle and number of photons measured by the detector of internally reflected 
Cherenkov light (DIRC). We accept K° s — ► 7r + 7r candidates that have a two-pion invariant mass 
within 12 MeV/c 2 of the nominal K® mass [13], a decay length greater than three standard devia- 
tions, and an angle between the line connecting the B and K® decay vertices and the direction of 
Kg momentum vector less than 45 mrad. The three daughters in the B decay are fitted constraining 
their paths to a common vertex, and the t: + tt~ mass to the nominal K® mass. 

In the characterization of the B candidates we use two kinematic variables. The energy differ- 
ence AE = Eb* — y/s/2 is reconstructed from the energy of the B candidate Eb* and the total en- 
ergy y/s, both evaluated in the e + e~ center-of-mass (CM) frame. The AE resolution for signal B de- 
cays is 18 MeV. We also use the beam-energy-substituted mass itt-es = J (s/2 + pi ■ pB) 2 /Ef — p|, 
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where (pi,Ei) is the four-momentum of the initial e + e system and ps is the momentum of the 
B candidate, both measured in the laboratory frame. The toes resolution for signal B decays is 
2.6 MeV/c 2 . We retain candidates with |A£| < 200 MeV and 5.2 < m E s < 5.3 GeV/c 2 . 

The background is dominated by random combinations of tracks created in e + e~ — ► qq (q = 
u, d, s, c) continuum events. We suppress this background by utilizing the difference in the topology 
in the CM frame between jet-like qq events and spherical signal events. The topology is described 
using the angle Qt between the thrust axis of the B candidate and the thrust axis of the charged 
and neutral particles in the rest of the event (ROE) [12]. Other quantities that characterize the 
event topology are two sums over the ROE: Lq = J2\Pi*\ an d £2 = J2 \P~i*\ cos 2 ®i-> where 6i is 
the angle between the momentum p* of particle % and the thrust axis of the B candidate. Addi- 
tional separation is achieved using the angle Ob between the S-momentum direction and the beam 
axis. After requiring \ cos 9t\ < 0.9, these four event-shape variables are combined into a Fisher 
discriminant T [14]. 

The remaining background originates from B decays where a neutral or charged pion is missed 
during reconstruction (peaking B background). We use samples of exclusive Monte Carlo (MC) 
events to model the signal and the peaking background, and data sidebands to model the continuum 
background. 

We suppress background from B decays that proceed through at^c transition leading to the 
K + K~K® final state by applying invariant mass cuts to remove D°, J/ip , Xc0> an d ip(2S) decaying 
into K + K~, and D + and Df decays into K + K®. Finally, to suppress B decays into final states 
with pions, we apply particle identification criteria to limit the pion misidentification rate to less 
than 2%. 



3 MEASUREMENT OF CP CONTENT 

The extraction of the CP content of the K + K~K® final state is based on an unbinned extended 
maximum likelihood fit. The likelihood function C is defined as: 



£ = exp II 



(2) 



where j runs over all 27368 K + K~K® events in the sample. The probability density function 
(PDF) V is formed as "P(toes) ■V(AE) - ViTY Event yields N\ for signal, continuum, and peaking 
B background are floated in the fit. We find a total of 525 ±30 signal events in the entire K + K~K® 
sample (including 4>K®), and show projection plots of the fit variables in Figure 1. 

The CP content is extracted using an angular moment analysis [15] which examines the distri- 
bution of the cosine of the helicity angle 6h between the K + and B° directions in the K + K~ center 
of mass frame. In this approach, we assume the decay rate for a given m(K + K~) invariant mass 
can be represented in terms of moments (Pi) of Legendre polynomials Pi(cos9h) (see Appendix A) 

\A\ 2 = ]r<P,).p,(cos<M, (3) 
1 

where A is the decay amplitude. Since the dynamics of the decay is not known, we extract the 
moments by summing over all events 

(Pt) « ^2 p i( cosd Hj) yvj/ej, (4) 
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Figure 1: Projection plots of the twes an d APJ variables. The points are data and the curves are 
projections from the likelihood fit for all events (full line) and continuum background (dashed line). 
The signal-to-background ratio is enhanced with a cut on the signal probability. 



where W is the weight for event j to belong to the signal decay. We use the s Plot background- 
substraction technique [16] to compute the weights as Wj = »7't>' ~ > wnere V 8 % is the signal row 
of the covariance matrix obtained from the fit. The efficiency e is evaluated from a high-statistics 
MC sample in m(KK) — cos Oh bins. The covariance matrix for the moments is computed as 

4 ~ £^(cos0tfj)fMcos0Hj) Wj/4 (5) 

j 

where the sum runs over events. Limiting ourselves to the two lowest partial waves, we can write 
the total decay amplitude in terms of the S-wave (CP-even) and the P-wave (CP-odd) amplitudes, 

A = A s P (cos6 H ) + e^ApP^cosOH), (6) 

where 4> p is the relative phase between the partial- wave amplitudes A s and A p . It can be easily 
shown that the angular moments (Po,2) give infomation on S- and P-wave strengths, and (Pi) arises 
from their interference. Comparing Equations (3) and (6), we can relate the moments with the 
wave intensities and the total fraction of CP-even events, f even as 



V2(P )-^(P 2 ), (7) 

k = yi^), (8) 



hven ~ Al + Al ~ V4(P ), (J) 

where A 2 and A p are the S- and P-wave intensities, respectively. As a cross-check, we extract the 
fraction of CP-even events by comparing the event rates of two isospin-equivalent channels [10]: 



sm T(B+ K+K° S K° S ) _ N K+K o K o {e K+K - K%/ Tj?0 

T(B° — > K + K~Kf) N K+K - K% (e K+KOgK o) r B+ { ] 

where efficiencies (e) are averaged over the Dalitz plot and include branching fractions for K® 
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4 MEASUREMENT OF CP ASYMMETRY 



The CP asymmetry parameters are extracted from a K + K~K® sample that excludes (f)K® decays 
by requiring \m(K + K~) — m{4>)\ > 15MeV/c 2 . We use the maximum-likelihood fit from Eq. (2), 
where the total PDF is formed as 'P(toes) ■ V(AE) • V(T) • V c (At; a At)- The time difference At is 
extracted from the measurement of the separation Az between the Bcp and £?tag vertices, along the 
boost axis (z) of the BB system. The vertex position of the Bcp meson is reconstructed primarily 
from the kaon tracks, and its MC-estimated resolution ranges between 40-80 /im, depending on the 
opening angle and direction of the kaon pair. The final resolution is dominated by the uncertainty 
on the #tag vertex which allows for a At (Az) precision with an r.m.s. of 1.1 ps (180 /xm). We 
retain events that have |Ai| < 20 ps and whose estimated uncertainty a At is less than 2.5 ps. The 
At resolution function is parameterized as a sum of two Gaussian distributions whose widths are 
given by a scale factor times the event-by-event uncertainty a At- A third Gaussian distribution, 
with a fixed large width, accounts for a small fraction of outlying events [17]. 

Recoil side decay products are used to determine the flavor of the -Stag m eson (flavor tag) and 
to classify the event into seven mutually exclusive tagging categories [17]. If the fraction of events in 
category c is e c and the mistag probability is w c , the overall quality of the tagging, J2 C e c(l — 2u> c ) 2 , 
is (30.5 ±0.5)%. Parameters describing the tagging performance and the At resolution function for 
signal events are extracted from approximately 30,000 B° decays into D ( -*^X + (X + = ir + , p + , af) 
flavor eigenstates. 

In the fit we float CP asymmetry parameters, parameters describing the At resolution function 
and tagging, event yields, and the signal PDF parameters for ttieSj AE, and the Fisher discriminant. 

5 SYSTEMATIC STUDIES 

In the measurement of the CP-even fraction based on the angular moments we estimate a bias due 
to the efficiency modeling from high-statistics MC events (2.5%). We do not find indication for 
the existence of higher moments (Pi), I = 3... 6 (see Figure 4 in Appendix A), that could arise 
from intermediate D-wave decays into K + K~ or decays proceeding through an / = 1 resonance 
into K^Kg. Nevertheless, we estimate a systematic error from the D-wave by examining the (P2) 
moment in the /2(1270), 0^(1320) and ^(1525) mass region (1.1-1.7 GeV/c 2 ) and assuming that 
(P2) arises only from D-wave and S-D interference. Since the moment itself is consistent with zero, 
we assign a systematic error of 4% based on the (P2) error. We also assign a systematic error 
due to potential contributions from decays proceeding through isovector resonances into K + K®. 
In addition to not being included in the angular moment analysis, decays of charged a^(980), 
a~2 (1320) and 0^(1540) are not CP eigenstates. We set the error on the CP content by counting 
events in K^K® invariant mass regions under these resonances [13]. We do not observe events 
in the ao(980) region, or a peak consistent with the ao(1540) resonance. We estimate events that 
could come from 02(1320) decays and conservatively include them into the systematic error (4.6%). 

The systematic errors on the time-dependent CP-asymmetry parameters (as, ac) are estimated 
similarly to our previous analysis [11]. We account for the fit bias (0.02, 0.01), the presence of 
double CKM-suppressed decays in -Btag [18] (0.018, 0.053), the uncertainty in the beam spot 
and detector alignment (0.022, 0.012), and the asymmetry in the tagging efficiency for signal 
and background events (0.011, 0.014). Other smaller effects come from the At resolution and 
uncertainty on the B° lifetime and mixing frequency (0.006, 0.006). We use t b o = 1.536 ±0.014 ps 
and Am d = 0.502 ± 0.007 ps^ 1 [13]. 
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Figure 2: Distributions of S- and P-wave intensities and CP even fraction as a function of K + K~ 
invariant mass. Insets show S- and P-wave intensities in the (f> mass region. Events within 15 MeV/c 2 
of the nominal (j) mass [13] are removed from the CP fit. 



6 RESULTS 

Distributions of the S- and P- wave intensities, and the CP-even fraction as a function of K + K~ 
invariant mass are shown in Figure 2. 4>K® events give a significant enhancement of P-wave decays 
in the first bin. In the sample that excludes c/)Kg events, we compute (Po,2) moments from the 
remaining events and find the total fraction of CP-even final states to be 

Uven = 0.89 ± 0.08 ± 0.06, 

where the first error is statistical and the second is systematic. We cross-check this result using 
the isospin approach of Eq (10). We find 452 ± 28 signal events in the K + K~ CP sample with 
a total efficiency of (e) = (17.3 ± 0.3)% and 208 ± 18 signal events in the K + K®K® sample with 
(e) = (9.8 ± 0.8)%. This gives feverf 1 = 0.75 ± 0.11 which is consistent with our nominal estimate. 

The coefficients of the time-dependent CP asymmetry in B° — > K + K~K® decays (excluding 
4>Kg final states) are determined to be 

S = -0.42 ±0.17 ±0.04, 
C = 0.10 ±0.14 ±0.06. 

The At distributions of events with B° and B° tags, with projections from the likelihood fit 
superimposed, are shown in Figure 3. The fit procedure is verified with the K + K®K® sample, 
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where we measure zero asymmetry, and the J/tpK® sample where the results are consistent with 
previous measurements [4, 5]. 




Figure 3: Distributions of At for B° — ► K + K~K® candidates with (a) B° tags and (b) B° tags. 
The solid lines refer to the fit for all events and the dashed lines correspond to the background 
contribution. The distribution of the raw asymmetry is shown in (c), where the solid line is obtained 
from the fit. The signal-to-background ratio is enhanced with a cut on the signal probability. 



The presence of both P- and S-wave decays in our CP sample dilutes the measurement of the 
sine coefficient. If we account for the measured CP-odd fraction, we can extract the SM parameter 
sin 2/3. Using the estimate of the CP content based on the angular moments, and setting C = in 
the fit, we get 

sin 2(5 = -S/{2f even - 1) = 0.55 ± 0.22 ± 0.04 ±0.11, 
where the last error is due to uncertainty on the CP content. 

7 SUMMARY 

In a sample of 227 million BB mesons, we have obtained preliminary measurements of the CP 
content and CP parameters in the K + K~K® final state that excludes (fiK® decays. From the 
distribution of the helicity angle in the K + K~ frame, described in terms of moments of Legendre 
polynomials, we extract the fraction of P-wave decays. The result is consistent with our cross- 
check and previous measurements based on isospin symmetry [7, 11], and confirms the dominance 
of CP-even final states. 

We measure a time-dependent CP asymmetry in B° — ► K + K~K® decays at the 2.3<r level. 
The obtained value for sin 2(3 is consistent with the SM expectation and previous measurements in 
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decays into the K + K K® final state [7, 11]. 
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which are orthogonal and normalized to unity J cI(cos9h)PiPi' = &w- Moments of Legendre poly- 
nomials (Pi) = J d(cos 9h)Pi(cos 9h)\A\ 2 can be extracted by replacing the integration over the 
unknown amplitude A with a sum over signal weights as shown in Eq. (4). The moments can be 
related with wave amplitudes A S;P;( i and interference phases 4>p,d as follows 



A APPENDIX: ANGULAR MOMENTS 



In the analysis we use the following definitions for the Legendre polynomials 
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(2 

di 



where we kept only terms relevant for the S, P and D waves (in the nominal result we set = 0). 
Lowest moments plotted in Figure 4 are used in the extraction of the P-wave fraction. Some of the 
higher moments shown in Figure 4 are used for systematic studies. In the normalization we assume 
that y/2 (Pq) equals the total number of signal events. 
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Figure 4: Distributions of angular moments (Pi) with I = 0, ... ,6. Insets show moments in the 
mass region. 
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